Phase Diagrams, Distinct Conformal Anomalies and Thermodynamics of Spin-1 Bond- Alternating 
Heisenberg Antiferromagnetic Chain in Magnetic Fields 
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The ground state phase diagrams and thermodynamic properties of a spin-1 bond-alternating Heisenberg anti- 
ferromagnetic chain with a single-ion anisotropy in longitudinal and transverse magnetic fields are investigated 
jointly by means of the infinite time evolving block decimation, the linearized tensor renormalization group 
and the density matrix renormalization group methods. It is found that in the magnetic field-bond alternating 
ratio plane six phases such as singlet dimer, Haldane, two Tomonaga-Luttinger liquid, 1/2 magnetization plateau 
and spin polarized phases are identified in a longitudinal field, while in a transverse field there are fives phases 
including singlet dimer, Haldane, Z 2 symmetry breaking, quasi- 1/2 magnetization plateau and quasi spin polar- 
ized phases. A reentrant behavior of the staggered magnetization in a transverse field is observed. The quantum 
critical behaviors in longitudinal and transverse magnetic fields are disclosed to fall into different universality 
classes with corresponding conformal field central charge c = 1 and 1 /2, respectively. The experimental data of 
the compound NTENP under both longitudinal and transverse magnetic fields are nicely fitted, and the Luttinger 
liquid behavior of low temperature specific heat experimentally observed is also confirmed. 
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I. INTRODUCTION 

Low dimensional quantum spin systems have been active 
subjects in quantum many-body physics for many years. Ow- 
ing to strong quantum fluctuations and competitions between 
various interactions in these systems, a number of exotic and 
fascinating quantum emergent phenomena are expected to oc- 
cur, which thus arouses persistently considerable interest not 
only in condense matter physics but also in other fields such 
as quantum information and quantum computation 1 . 

Among others, the spin-1 bond alternating Heisenberg an- 
tiferromagnetic chain (BAHAFC) in longitudinal and trans- 
verse magnetic fields is of particular interest, for a series of 
materials, such as NENP 2 , NMOAP 3 , NDOAP 4 , and NTENP 5 
can be well described by the spin S=l BAHAFC, and the 
ideal model compound NiCC^^XNC^ClCU (NTENP) 
has been extensively studied both experimentally and theoret- 
ically in recent years [3.4.6-12], where the effects of magnetic 
field and the bond alternating ratio on low-lying magnetic ex- 
citations, spin correlations and low-temperature specific heat 
of the compound NTENP have been investigated. It is known 
that in a magnetic field the one-dimensional (ID) Heisenberg 
quantum spin chains with periodic ground states would ex- 
hibit a topological quantization of magnetization according to 
the necessary condition n(S - m) = integer—, where n is the 
period of the system, S is the spin and m is the magnetization 
per site. Such magnetic plateau states have been addressed 
in some polymerized Heisenberg antiferromagnetic or ferri- 
magnetic spin chains (e.g. Refs. 11411 ). High-field magnetiza- 
tion measurements on NTENP also revealed that am = 1/2 
magnetization plateau appears around 700 kOe 3,4 ; the neu- 
tron scattering experiments indicated that the spin dynamics 
of NTENP quite differs from that of a Haldane spin chai n 7 ' 8,12 ; 
and the low-temperature specific heat measurement showed 
a Tomonaga-Luttinger liquid (TLL) behavior in NTENP in 
a longitudinal magnetic field 11 . Although the experimental 
studies on NTENP gain obvious advances, some ambiguities 



still remain on further understanding its physical properties. 
For instance, the compound NTENP was experimentally stud- 
ied under both longitudinal and transverse magnetic fields, but 
only the data in a longitudinal field have been fitted by quan- 
tum Monte Carlo simulations; no direct numerical evidence of 
low-temperature specific heat of this material was reported to 
verify the TLL behavior; the sharp peaks of the specific heat 
in the NTENP were observed experimentally in both longi- 
tudinal and transverse magnetic fields, where the peak posi- 
tions move to high temperature side with the increase of the 
transverse field, while retain almost intact with increasing the 
longitudinal field; and so on. The remaining issues are very 
worthy to address. 

Apart from that the spin-1 BAHAFC with a single-ion 
anisotropy [Eq. ([TJ below] is believed to be a pertinent model 
in describing the physical characters of NTENP, which has 
been studied previouslyisliii, this model itself has also fasci- 
nating properties that deserve to investigate. When the bond 
alternating ratio a = 0, it reduces to the singlet dimers; when 
a - 1, it becomes a spin-1 uniform chain (Haldane chain). In 
these two special cases, the system shows an excitation gap 
from the singlet dimer or Haldane ground state to the triplet 
excited state. The latter uniform chain under both longitudinal 
(h z ) and transverse (h x ) magnetic fields has been extensively 
discussed (e.g. Refs. i2Tj| - |22ll ). However, the overall phase 
diagrams in a - h z x plane for the S=l BAHAFC with a single- 
ion anisotropy are still absent. Therefore, it is really necessary 
to tackle these above questions. 

In this paper, by means of jointly the infinite time evolv- 
ing block decimation (iTEBD) 23 , the linearized tensor renor- 
malization group (LTRG) 24 , and the density matric renormal- 
ization group (DMRG) methods^ 5 -, we shall study systemati- 
cally the ground state phase diagrams, magnetic and thermo- 
dynamic properties of the S=l BAHAFC with a single-ion 
anisotropy in longitudinal and transverse magnetic fields. The 
iTEBD, LTRG and DMRG methods are the numerical algo- 
rithms with very good accuracy and high efficiency that were 
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recently developed for low-dimensional quantum lattice sys- 
tems, which allow for calculating with nice precision the criti- 
cal properties and the extremely low temperature behaviors of 
quantum many-body lattice systems. These methods can as- 
sist us to fit well the experimental data on the NTENP with the 
accurately calculated results, thereby capable of reasonably 
determining the material parameters and better understanding 
the fundamental features of the NTENP in both longitudinal 
and transverse magnetic fields. 

II. MODEL AND GROUND STATE PHASE DIAGRAMS 
A. Model Hamiltonian and Calculational Method 

Let us start with the model Hamiltonian given by 

L/2 

H = /2j(S 2i -i • S 2i + aS 2j • S 2 ,+i) 

i 

L 

+ J][A(SZ) 2 -g ]} Li B h z SZ-g ± Li B h x Sfl (1) 

i 

where J is the coupling constant, L (even) is the number of lat- 
tice sites, S,- is the spin operator at /th site, A is the single-ion 
anisotropy, g\\, g ± are the Lande g-factors, /i B is the Bohr mag- 
neton, h z and h x are the longitudinal and transverse magnetic 
fields, respectively. 

In the following calculations, gyyUg = g±HB = 1 is assumed 
unless the fittings to experimental data are concerned, and J 
is taken as the energy scale. The ground state properties were 
studied by utilizing the iTEBD 23 imaginary time projection 
scheme along with DMRG method 25 , and the thermodynamic 
properties and fittings to the experimental data of NTENP 
were performed by invoking the LTRG approach^!. In iTEBD 
calculations we take the smallest Trotter step r = 1(T 7 . 

B. Magnetization and Phase Diagram in a Longitudinal Field 

Figure 1(a) presents the field dependence of magnetization 
per site m z for the system defined by Eq. (0 in the longitu- 
dinal magnetic field for a - 0.45 and 0.8, where the bond 
dimension D c = 120 of a matrix product state is presumed. 
When a = 0.45, it is clear to see that there are three mag- 
netization plateaux occurring at m=0, 1/2 and 1 in some re- 
gions of h z . The m=0 plateau suggests that a gap remains in 
the absence of magnetic field, implying that the system is in 
the singlet-dimer phase, as it should be in the same phase as 
a — 0. Based on the valence-bond solid picture, the m =1/2 
plateau corresponds to that one of the two bonds in each cou- 
pled spin-1 dimer is broken 14 . The m = 1 plateau is in the spin 
fully polarized state. Near the critical fields, the magnetiza- 
tion per site m z depends on the magnetic field in a square root 
behavior, similar to the case of spin-1 HAFC in a magnetic 
field. In the ranges between the plateaus, the spin correlations 
should be in a power law decay, indicating a quasi long-range 
order. The critical magnetic fields for closure of the gap can 
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FIG. 1: (Color online) (a) The magnetic curves of spin-1 bond alter- 
nating Heisenberg antiferromagnetic chain (BAHAFC) with a single- 
ion anisotropy in longitudinal magnetic fields h z , where three mag- 
netization plateaux (m = 0, 1/2, and 1) are observed. The inset 
shows the singularities of the susceptibility as a function of magnetic 
field, which signal the occurrence of quantum phase transitions, (b) 
The ground state phase diagram in the plane of the bond alternating 
ratio (a) vs. longitudinal magnetic field (h z ), where the Haldane, 
singlet-dimer (SD), m=l/2 magnetic plateau, Tomonaga-Luttinger 
liquid (TLL), and spin polarized phases are identified. The single- 
ion anisotropy A = 0.25 is taken. 



be accurately determined by the field dependence of the sus- 
ceptibility, where the sharp peaks appear, as shown in the inset 
of Fig. 1(a). When a = 0.8, the three magnetic plateaux still 
remain, but the width of 1/2 plateau decreases, indicating that 
the corresponding gap becomes narrower. For the spin-1 BA- 
HAFC, as the ground state period is n = 2, the topological 
quantization condition n(S - m) = integer is satisfied by m=0, 
1/2 and 1, i.e., there should exist at most three plateaux in the 
longitudinal field for < a < 1, which is well confirmed in 
Fig. 1(a). At the special case of a — 1, n becomes 1, and there 
should be two plateaus at m-Q and 1 . 

ForO < a < 1, we have calculated the magnetization curves 
and the corresponding susceptibilities in the ground state, and 
collected all critical values of magnetic fields at which the sus- 
ceptibility exhibits obvious singularities, which allows us for 
readily obtaining a global ground state phase diagram in a— h z 
plane for the spin-1 BAHAFC with a single-ion anisotropy 
A = 0.25, as shown in Fig. 1(b). There are six phases, includ- 
ing the singlet-dimer (SD) phase, Haldane phase, two TLL 
phases separated by the m = 1/2 magnetic plateau phase, 
and the spin fully polarized phase. At a critical bond alter- 
nating ratio a c ~ 0.6, there exists a quantum phase transition 
from the gapped SD state into the gapped Haldane phase. The 
m-l/2 plateau state persists for all values of a until it disap- 
pears at a = 1. At all six phase boundaries, the gap closes. 
In these phases, the spin-spin correlation function (S^S^) re- 
veals different spatial behaviors, which is short-range ordered 
except that in the TLL phase it decays in an algebraic way. 
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FIG. 2: (Color online) (a) The magnetic curves of S = l BAHAFC 
with a single-ion anisotropy in transverse magnetic fields h x . The up- 
left inset shows the staggered magnetization along the y axis, and the 
down-right inset presents the corresponding susceptibility as a func- 
tion of magnetic field. For a = 0.1, three plateau-like steps are ob- 
served, (b) The ground state phase diagram in the plane of the bond 
alternating ratio (a) vs. transverse magnetic field (h x ), where the Hal- 
dane, singlet-dimer (SD), Z2 symmetry breaking Neel ordered, quasi 
1/2 magnetic plateau, and quasi spin polarized phases are identified. 
The single-ion anisotropy A = 0.25 is taken. 



C. Magnetization and Phase Diagram in a Transverse Field 

In a transverse magnetic field h x , the field dependence of 
magnetization per site m x and m y and the corresponding sus- 
ceptibility of the model under interest are presented in Fig. 
2(a) for a — 0.1 and 0.45. Three steps around m x =0, 1/2 and 
1 in magnetic curves (m x vs. h x ) are observed for a = 0.1. 
The corresponding susceptibility x is illustrated in the lower 
inset of Fig. 2(a), where one may notice that, due to the non- 
conservation of total S x , at these three steps \ does not exactly 
vanish although it is negligibly small. Thus, we call these 
steps as quasi magnetic plateaus. In the regions between the 
quasi plateaus, m y has large values as shown in the upper inset 
of Fig. 1(a). The reason is that the existence of the single-ion 
anisotropy makes the spins tend to arrange in xy plane, and 
when the magnetic field is applied along the x direction, the 
staggered magnetization m y is induced to form a canted Ising 
order. In addition, for a = 0.1 we observed from the upper 
inset of Fig. 2(a) that with increasing h x the staggered magne- 
tization m y first is negligible small, and then exhibits a sharp 
valley structure, which illustrates a reentrant behavior of the 
staggered magnetization 26,27 . For a = 0.45, only two steps 
appear in the magnetic curve, and the quasi m = 1/2 step dis- 
appears. In this case, m y has even larger values in the region 
between the two steps. For other larger a, the magnetic curves 
in the transverse field h x have the behaviors similar to that of 
a = 0.45. 

Utilizing the method similar to Fig. 1(b) and sweeping var- 
ious values of a, we can collect all critical magnetic fields by 
finding the singular positions of susceptibility, and then obtain 
the whole ground state phase diagram in a - h x plane for the 
system under investigation, as depicted in Fig. 2(b). It can 
be seen that there are five phases, namely, the SD, Haldane, 
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FIG. 3: (Color online) The normalized Schmidt coefficient A v as a 
function of the number of kept states v for different phases, where 
D c = 80. 



quasi m = 1/2 magnetic plateau, Z2 symmetry breaking (or 
spin canted) and the quasi spin polarized phases. The quan- 
tum critical point a c =s 0.6 is recovered here. In contrast to 
the case under a longitudinal field where the m = 1/2 plateau 
phase persists into the whole region of a < 1 and separates 
two TLL phases, the quasi m = 1/2 plateau phase sets only in 
a small region (a < 0.24). Starting from the SD and Haldane 
phases, when we increase the field h x , the gap is gradually 
closed at the phase boundaries. If one continues to increase 
the magnetic field h x , another gap opens because of the break- 
ing of the discrete symmetry, giving rise to a Z2 symmetry 
breaking phase, which displays a long range behavior. 



III. ENTANGLEMENT ENTROPY AND CONFORMAL 
ANOMALIES 



This present system demonstrates quite different behaviors 
in longitudinal and transverse magnetic fields. To gain fur- 
ther insight into the underlying physics behind this character, 
we have studied the conformal anomalies of this system at the 
critical regimes. The conformal field theory (CFT) tells us 
that the conformal invariance at the critical point sets useful 
constraints on the critical behaviors of two-dimensional clas- 
sical or ID quantum systems 28 , and the universality class can 
be characterized by the conformal anomaly or central charge c 
of the Virasoro algebra. For a system with a continuous sym- 
metry G, if G=SU(N), the possible conformal central charge 
is given by c = (N 2 - l)k/(N + k), k = 1, 2, 3, ..M. For the 
spin-S uniform antiferromagnetic quantum chains, G=SU(2), 
c = 35/(1 + S). For S=l/2, c = 1; S=l, c = 3/2. For 
the present spin-1 BAHAFC with a single-ion anisotropy, the 
SU(2) symmetry is not satisfied, and the conformal central 
charge should be calculated via other ways. The von Neu- 
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mann entropy S , defined by 

S = -Tr{p sys \og 2 p sys ) = -Tr(p mv log 2 p env ), (2) 

offers a possible way to get the central charge of quantum spin 
chains, where p sys ( e nv) is the reduced density matrix (DM) of 
system (environment). In critical and noncritical regimes, the 
entanglement entropy has different asymptotic behavior a 31 ' 32 . 
In critical regimes, the CFT predicts^ 3 - 

S *|log 2 (L)+fc, (3) 

where L is the number of spins for a block embedding in an 
infinite chain, c is the central charge and A: is a non-universal 
constant. In noncritical regimes, S vanishes for all L or grows 
monotonically with L until saturation. In the following, we 
shall invoke the iTEBD method to study the entanglement en- 
tropy for the case in a transverse magnetic field. Because of 
the nearly unitary evolution, the iTEBD algorithm gives the 
canonical form of infinite matrix product states (iMPS). We 
can make a Schmidt decomposition on the ground state wave 
function such that 

m = 2^AAtf)\tf) (4) 

v=l 

on each bond, where \if/y' B ) is the orthonormalized basis states 
(Schmidt bases) for two parts A and B of the infinite chain, and 
A v is the Schmidt coefficient (SC). Fig. [3]presents A v as a func- 
tion of v for different phases. It is seen that in the TLL phase, 
the SC attenuates much slower than those in other gapped 
phases. The double degeneracy of SC only appears in the Hal- 
dane phase, indicating the existence of the topological string 
order—. Given a canonical form of iMPS, if the system in Eq. 
(O is chosen as the semi-infinite chain, it is easy to directly 
calculate the von Neumann entropy, S - - log 2 (/lJ). 
In addition, we can also calculate the entanglement entropy 
when the system is successive L spins embedding in an in- 
finite chain, where we may set L = 2 for an example. In 
particular, we can write 

D c d 

W>= Yj Z ^Kfi^iy^MjWr), (5) 

aJ3,y=l i,]=\ 

where P aj3 is the tensor on the i-th site, is the spin basis 

states on the /(y')-th site, \<i) a ) and |O r ) are the left Schmidt 
basis of the i-th site and the right Schmidt basis of the j-th 
site, respectively. By contracting two tensors, we will obtain 
the reduced DM 

p = TryMffl 

D c D c d D c 

= Zj Z [ Z Z AaV «fi A P Y iy A y 

a,y=\a',y'=\ i,j=l f3J3'=l 

^(rl^%'(XyyrAy]\a)\y)(a'\(y'l (6) 

One can diagonalize this matrix and calculate the entangle- 
ment entropy. For a larger L, we need to contract L tensors 
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FIG. 4: (Color online) The entanglement entropy 5 as a function 
of chain length L calculated by the iTEBD in a transverse field for 
a = 0.45 and D c = 100, and by the DMRG in a longitudinal mag- 
netic field for a = 0.6 where the optimal states were kept as 400. In 
both cases, the single-ion anisotropy A = 0.25. The solid line is the 
fitting curve to Eq. (3), giving the central charge c = 1/2 at a critical 
transverse field h\ = 0.349; the inset shows that the entanglement 
entropy as a function of h x for a semi-infinite chain length is singular 
at h c x = 0.349. The dashed line is the fitting curve for the case (TLL 
phase) at a longitudinal magnetic field h- = 0.3, giving the central 
charge c = 1 . 



and select Schmidt basis to construct the reduced DM. For all 
L, the dimension of the reduced DM is D 2 C x D 2 C . By calcu- 
lating the von Neumann entropy for a semi-infinite chain in a 
transverse magnetic field, we observe that the cusp position of 
the entanglement entropy just gives the phase transition point 
h\ = 0.349(2). By fitting our calculated results to Eq. d3), we 
find c = 1/2 for h% = 0.349, as shown in Fig. H In the TLL 
phase of the case in a longitudinal magnetic field, the iTEBD 
algorithm gives the ground state with a staggered magneti- 
zation perpendicular to z direction so that it cannot yield a 
correct wave function. The reason is that there does not have 
an excitation gap and the iTEBD algorithm cannot project the 
right ground state wave function. In order to calculate the 
central charge in the TLL phase, we utilize the DMRG algo- 
rithm with open boundary condition at h z = 0.3, a - 0.6, 
and A = 0.25. The result is also included in Fig. [4] where 
the fitting result gives the central charge c = 1 in this TLL 
phase 35 . Consequently, this implies that the critical behaviors 
of the spin-1 BAHAFC with a single-ion anisotropy in trans- 
verse fields quite differ from those in longitudinal fields, and 
the universality falls into two distinct classes with conformal 
central charges c = 1/2 and c = 1, respectively. 
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FIG. 5: (Color online) The temperature dependence of susceptibil- 
ity of NTENP measured experimentally is well fitted to the LTRG 
calculated data for both longitudinal and transverse magnetic fields. 
The experimental data are taken from Ref. |4j|. 



IV. THERMODYNAMIC PROPERTIES AND 
COMPARISON TO EXPERIMENTS 

A. Susceptibility, Magnetization and Comparison to 
Experiments 

The thermodynamic properties of this system were ex- 
plored by using recently proposed LTRG algorithm 24 that al- 
lows for accurately calculating the thermodynamic quantities 
at very low temperature. In the following LTRG calculations 
we keep the Trotter step r = 0.1. 

Fig. [5] gives the temperature dependence of susceptibil- 
ities of the spin-1 BAHAFC in longitudinal and transverse 
magnetic fields, where the fittings to the experimental data 
of NTENP under both fields are also included. One may see 
that the theoretical results are nicely fitted to the experimen- 
tal data, generating a set of material parameters of NTENP: 
a = 0.45, A/7 = 0.25, J = 54.2K, g {{ = 2.14, and g ± = 2.27. 
It should be noted that in a longitudinal field, these parameters 
are in agreement with the previous results from the quantum 
Monte Carlo calculations 4 , while in a transverse field the fit- 
ting to the experimental data of NTENP is for the first time 
done. The magnetization curves of NTENP were measured 
at T= 1.3 K up to 700 kOe in Ref. 0], which are well fitted 
to our LTRG calculated results, giving the same set of fitting 
parameters as those obtained from the data of susceptibility 
except that g ± = 2.24 here, as shown in Fig. [6] Consider- 
ing that the two sets of experimental data from susceptibility 
and magnetic curves are independent, a slight difference on 
g ± is reasonable. Thus, g ± of NTENP should be around 2.24 
~ 2.27. The high-field magnetization curve is also well fit- 
ted with our LTRG results (inset of Fig. ((6}) using the same 
parameters. 
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FIG. 6: (Color online) The magnetic curves of NTENP are well fit- 
ted to the LTRG calculated data for both longitudinal and transverse 
fields, where the fitting parameters are consistent with those from the 
susceptibilities. Inset is the high-field magnetic curve up to 700 kOe 
fitted with the corresponding LTRG calculated data. The experimen- 
tal data are taken from Ref. |4fl. 



B. Specific Heat in Magnetic Fields 

The temperature dependence of the specific heat C(T) of 
the S=l BAHAFC in longitudinal and transverse magnetic 
fields is obtained by the LTRG method down to very low tem- 
peratures, as shown in Fig. Qfor h z = 0.5 and h x = 0.5. In 
a transverse field, there are one round peak and a low tem- 
perature shoulder in C(T), while in a longitudinal field the 
specific heat exhibits only one broad peak. It appears that 
at low temperature the specific heat displays distinct behav- 
iors. To show this point clearly, we have carefully calculated 
the specific heat of this model at extremely low temperatures 
by the LTRG algorithm that is very powerful and efficient for 
calculating the low-temperature properties than other numer- 
ical methods 36 . Shown in the lower inset of Fig. [7] two dif- 
ferent behaviors in both fields at low temperature are clearly 
demonstrated, where in a longitudinal field, the specific heat 
displays a linear T-dependent character, showing a TLL be- 
havior, while in a transverse field, C(T) exhibits an exponen- 
tial decay (the fitting curve is given in the upper inset of Fig. 
|7J, that can be ascribed to the Z2 symmetry breaking with an 
open of an Ising gap. 

The specific heat C(T) of NTENP was also measured 
experimentally-ii, where the linear T-dependence of C(T) at 
low temperature in a longitudinal field above the critical field 
h z ° = 9.3 T was observed, which was identified as a TLL be- 
havior. Our low temperature LTRG calculations at h z = 0.5 
(> h z c ) on the S=l BAHAFC model strongly supports this 
experimental observation. In a transverse field, the experi- 
ment on NTENP gives a distinct nonlinear low-temperature 
behavior of C(T) from that in a longitudinal field, which is 
also backed up by our LTRG results. We should remark here 
that the sharp peaks in the temperature dependent specific 
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FIG. 7: (Color online) The temperature dependence of the specific 
heat of S=l BAH AFC under both longitudinal and transverse mag- 
netic fields. Although the whole profiles look similar for both cases, 
the extremely low temperature behaviors shown in the down-left in- 
set quite differ, where the linear-T dependence of the specific heat in 
a longitudinal field is clearly seen at low temperature, suggesting a 
Tomonaga-Luttinger liquid behavior, while an exponential decay for 
the T-dependence of the specific heat is observed in a transverse field 
owing to the appearance of an Ising gap, as seen in the up-right inset. 
The single-ion anisotropy A = 0.25 and the bond alternating ratio 
a = 0.45 are taken. 



heat of NTENP were observed in both fields, which cannot 
be explained by using this spin-1 BAHAFC model, as re- 
vealed by our present studies, because we do not find any 
field-dependent sharp peaks of C(T) in this ID model. Those 
sharp peaks may signal the field-induced long-range orders 
from the 3D effect of inter-chain interactions in NTENP. It 
is this reason that makes us not directly fit our LTRG results 
with the experimental data of specific heat of NTENP. Nev- 
ertheless, our present studies on the spin-1 BAHAFC model 
may give a possible clue to understand the experimentally ob- 
served low-temperature sharp peaks of C(T)^. In a longitudi- 
nal field, as long as the magnetic field is higher than the criti- 



cal field (h z c ), the system will enter the TLL phase that has a 
quasi LRO, and the true LRO can be established only by inter- 
chain couplings. Thus the peak position is almost determined 
by the strength of interchain interactions and hardly moves 
with the increase of the magnetic field. In a transverse field, 
in the range that the experiment was performed, the staggered 
magnetization that could enhance the interchain couplings in- 
creases monotonously with increasing the magnetic field, so 
the low-temperature sharp peak moves to the high tempera- 
ture side with the increase of the transverse field. 



V. SUMMARY AND CONCLUSION 

In summary, we have investigated the spin-1 BAHAFC 
model with a single-ion anisotropy in longitudinal and trans- 
verse magnetic fields by employing the iTEBD and LTRG 
methods. The ground state phase diagrams in the plane of 
field versus the bond alternating ratio under both fields are 
obtained, where various phases are identified for two cases. 
From the entanglement entropy, the conformal central charges 
in both critical longitudinal and transverse fields are deter- 
mined to be c = 1 and 1 /2, respectively, suggesting that the 
universality in critical regimes falls into different classes for 
both fields. The TLL behavior at low T observed experimen- 
tally in NTENP is verified via our accurate calculations, while 
an exponential decay of low-temperature specific heat is un- 
covered in a transverse field. The experimental data of the 
model material NTENP are well fitted with our LTRG results, 
and the parameters for characterizing NTENP are determined . 



Acknowledgments 

We are indebted to Bo Gu, Fei Ye, Shou-Shu Gong, Bin 
Xi, J. Sirker, and Qing-Rong Zheng for stimulating discus- 
sions. This work is supported in part by the NSFC (Grants 
No. 10934008, No. 90922033), the MOST (Grant No. 
2012CB932900) and the Chinese Academy of Sciences. 



* Corresponding author. Email: gsu@gucas.ac.cn 
* These authors contribute equally to this work. 

1 A. S. Darmawan, and S. D. Bartlett, Phys. Rev. A 82, 012328 
(2010). 

2 S. Ma, C. Broholm, D. H. Reich, B. J. Sternlieb, and R. W. Erwin, 
Phys. Rev. Lett. 69, 3571 (1992) 

3 Y. Narumi, M. Hagiwara, R. Sato, K. Kindo, H. Nakano, and M. 
Takahashi, Physica B 509, 246-247 (1998). 

4 Y. Narumi, K. Kindo, M. Hagiwara, H. Nakano, A. Kawaguchi, 
K. Okunishi, M. Kohno, Phys. Rev. B 69, 174405 (2004). 

5 A. Escuer, R. Vicente, and X. Solans, J. Chem. Soc. Dalton Trans. 
(1997), 531. 

6 Y. Narumi, M. Hagiwara, M. Kohno, and K. Kindo, Phys. Rev. 
Lett. 86, 324 (2001). 

7 A. Zheludev, T. Masuda, B. Sales, D. Mandrus, T. Papenbrock, T. 
Barnes, and S. Park, Phys. Rev. B 69, 144417 (2004). 



8 M. Hagiwara, L. P. Regnault, A. Zheludev, A. Stunault, N. 
Metoki, T. Suzuki, S. Suga, K. Kakurai, Y. Koike, P. Vorderwisch, 
and J. -H. Chung, Phys. Rev. Lett. 94, 177202 (2005). 

9 T. Suzuki and S. Suga, Phys. Rev. B 72, 014434 (2005). 

10 L. P. Regnault, A. Zheludev, M. Hagiwara, and A. Stunault, Phys. 
Rev. B 73, 174431 (2006). 

11 M. Hagiwara, H. Tsujii, C.R. Rotundu, B. Andraka, Y. Takano, N. 
Tateiwa, T.C. Kobayashi, T. Suzuki, and S. Suga, Phys. Rev. Lett. 
96, 147203 (2006). 

12 V. N. Glazkov, A. I. Smimov, A. Zheludev and B. C Sales, Phys. 
Rev. B 82, 184406 (2010). 

13 M. Oshikawa, M. Yamanaka, and I. Affleck, Phys. Rev. Lett. 78, 
1984(1997). 

14 B. Gu, G. Su, and S. Gao, J. Phys.: Condens. Matter 17, 6081 
(2005); B. Gu, G. Su, and S. Gao, Phys. Rev. B 73, 134427 (2006); 
B. Gu and G. Su, Phys. Rev. Lett. 97, 089701 (2006); Phys. Rev. 



7 



B 75, 174437 (2007); S. S. Gong and G. Su, ibid. 78, 104416 
(2008); S. S. Gong, S. Gao, and G. Su, ibid. 80, 014413 (2009); 
S. S. Gong and G. Su, Phys. Rev. A 80, 012323 (2009); Y. Zhao, 
S. S. Gong, W. Li, and G. Su, Appl. Phys. Lett. 96, 162503 (2010); 
S. S. Gong, W. Li, Y. Zhao, and G. Su, Phys. Rev. B 81, 214431 
(2010). 

15 G. Gomez-Santos, Phys. Rev. Lett. 63, 790 (1989). 

16 A. M. Tsvelik, Phys. Rev. B 42, 10499 (1990). 

17 R. R. P. Singh and M. P. Gelfand, Phys. Rev. Lett. 61, 2133 (1988). 

18 M. Kohno and M. Takahashi, Phys. Rev. B 57, 1046 (1998). 

19 A. Kitazawa, K. Nomura, and K. Okamoto, Phys. Rev. Lett. 76, 
4038 (1996). 

20 I. Affleck, Phys. Rev. B 43, 3215 (1991). 

21 Z. C. Gu and X. G. Wen, Phys. Rev. B 80, 155131 (2009). 

22 H. Z. Xing, G. Su, S. Gao, and J. H. Chu, Phys. Rev. B 66, 054419 
(2002). 

23 G. Vidal, Phys. Rev. Lett. 98, 070201 (2007); R. Orus and G. 
Vidal, Phys. Rev. B 78, 155117 (2008). 

24 W. Li, S.-J. Ran, S.-S. Gong, Y. Zhao, B. Xi, F. Ye, and G. Su, 
Phys. Rev. Lett. 106, 127202 (2011). 

25 S. R. White, Phys, Rev. Lett. 96, 2863 (1992). 

26 P. Azaria, H. T. Diep, and H. Giacomini, Phys. Rev. Lett. 59, 1629 
(1987). 



27 Y. Hieida, K. Okunishi, and Y. Akutsu, Phys. Rev. B 64, 224422 
(2001). 

28 A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov, J. Stat. 
Phys. 34, 763 (1984); Nucl. Phys. B 241, 333 (1984); D. Friedan, 
Z, Qiu, and S. Shenker, Phys. Rev. Lett. 52, 1575 (1984). 

29 I. Affleck, Phys. Rev. Lett. 56, 746 (1986). 

30 H. W. G. Blote, J. L. Cardy, and M. P. Nightingale, Phys. Rev. 
Lett. 56, 742(1986). 

31 G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. Lett. 90, 
227902 (2003). 

32 J. Eisert, M. Cramer, and M. B.Plenil, Rev. Mod. Phys. 82, 277 
(2010). 

33 C. Holzhey, F. Larsen, and F. Wilczek, Nucl. Phys. B 424, 443 
(1994). 

34 F. Pollmann, A. M. Turner, E. Berg, M. Oshikawa, Phys. Rev. B 
81, 064439 (2010). 

35 In the infinite chain DMRG algorithm with open boundary con- 
dition, the entanglement entropy in critical regimes should take 
the form of S L ~ (c/6) log 2 (L) + k'\ see also P. Calabrese, and J. 
Cardy, J. Stat. Mech. P06002 (2004). 

36 B. Gu, Ph. D. Thesis, Graduate University of Chinese Academy 
of Sciences, Beijing, China, (2006). 



